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ABSTRACT 

A complete, flux density limited sample of 96 faint (> 0.5 mJy) radio sources is selected 
from the IOC survey at 15.7 GHz in the Lockman Hole. We have matched this sample to a 
range of multi-wavelength catalogues, including SERVS, SWIRE, UKIDSS and optical data; 
multi-wavelength counterparts are found for 80 of the 96 sources and spectroscopic redshifts 
are available for 24 sources. Photometric reshifts are estimated for the sources with multi¬ 
wavelength data available; the median redshift of the sample is 0.91 with an interquartile range 
of 0.84. Radio-to-optical ratios show that at least 94 per cent of the sample are radio loud, 
indicating that the IOC sample is dominated by radio galaxies. This is in contrast to samples 
selected at lower frequencies, where radio-quiet AGN and starforming galaxies are present 
in significant numbers at these flux density levels. All six radio-quiet sources have rising 
radio spectra, suggesting that they are dominated by AGN emission. These results confirm 
the conclusions of Paper I that the faint, flat-spectrum sources which are found to dominate 
the IOC sample below ~ 1 mJy are the cores of radio galaxies. The properties of the IOC 
sample are compared to the SKADS Simulated Skies; a population of low-redshift starforming 
galaxies predicted by the simulation is not found in the observed sample. 
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1 INTRODUCTION 

Most studies of the faint (< 1 mJy) radio sky have focused on fre¬ 
quencies around 1.4 GHz, due to the increased telescope time re¬ 
quired to survey a field to an equivalent depth at a higher frequency. 
This means that while the composition of the lower-frequency ra¬ 
dio sky is well constrained (e.g. Padovani et al. 2009; de Zotti et al. 
2010; McAlpine et al. 2013), the faint radio sky at higher frequen¬ 
cies is relatively unstudied. The Tenth Cambridge (IOC; AMI Con¬ 
sortium: Davies et al. 2011 ; AMI Consortium: Franzen et al. 2011) 
survey at 15.7 GHz has covered x: 27 deg^ in ten different fields 
to a completeness limit of 1 mJy and a further 12 deg^ to a com¬ 
pleteness limit of 0.5 mJy. The IOC survey therefore provides the 
ideal starting point from which to study the faint, high-frequency 
sky. This survey has recently been extended to even fainter flux 
densities in two of the IOC fields by Whittam et al. (in prep), who 
calculated the 15 GHz source counts down to 0.1 mJy. 

To investigate the nature of the extragalactic radio source pop¬ 
ulation multi-wavelength studies are required, as the power-law 
nature of radio spectra means that radio data alone is not sufficient 
to classify source types or estimate redshifts. In the first paper in 
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this series studying the faint radio source population at 15.7 GHz 
(Whittam et al. 2013, hereafter Paper I) we discuss a sample of 
sources selected from the Lockman Hole for which data over a 
wide range of frequencies are available. The Lockman Hole is a 
region of the sky centred near 10’'45'", -1-58° (J2000 coordinates, 
which are used throughout this work) with exceptionally low Hi 
column density (Lockman et al. 1986). The low infrared back¬ 
ground (0.38 MJy sr“' at 100 pm; Lonsdale et al. 2003) in this 
area of the sky makes it ideal for deep extragalactic infrared ob¬ 
servations. As a result, as part of the Spitzer Wide-area Infrared 
Extragalactic survey (SWIRE; Lonsdale et al. 2003) sensitive in¬ 
frared observations of x: 14 deg^ of the Lockman Hole area have 
been made. The availability of deep infrared observations in the 
Lockman Hole has triggered deep observing campaigns at optical. 
X-ray and radio wavelengths (e.g. Biggs & Ivison 2006; Brunner et 
al. 2008; Ishisaki et al. 2001; Lonsdale et al. 2003; Mauduit et al. 
2012; Wright et al. 2010). The availability of data at such a wide 
range of frequencies makes the Lockman Hole a particularly good 
area for a multi-wavelength study of the faint radio source popula¬ 
tion. 

In Paper I we selected a sample of 296 sources from the IOC 
sample in the Lockman Hole. By matching this catalogue to sev¬ 
eral lower frequency surveys we have investigated the radio prop- 
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erties of the sources in this sample; all but 30 of the IOC sources 
are matched to sources in one or more of these surveys. We found 
a significant increase in the proportion of flat-spectrum sources at 
flux densities below s;! mjy - the median spectral index between 
15.7 GHz and 610 MHz changes from a = 0.75 for flux densit¬ 
ies greater than 1.5 mJy to a = 0.08 for flux densities less than 
0.8 mJy (the convention 5 cv v““, for a source with flux density 5 
at frequency v, is used throughout this work). This suggests that a 
population of faint, flat-spectrum sources is emerging at flux dens¬ 
ities < 1 mJy in the high-frequency sky. 

In Paper I the spectral index distribution of this sample of 
sources selected at 15.7 GHz was compared to those of two samples 
selected at 1.4 GHz from FIRST and NVSS. This showed that 
there is a significant flat-spectrum population present in the IOC 
sample which is missing from the samples selected at 1.4 GHz. 
The IOC sample was compared to a sample of sources selected 
from the Square Kilometre Array Design Studies (SKADS) Sim¬ 
ulated Sky (S^) by Wilman et al. (2008, 2010) and we found that 
this simulation fails to reproduce the observed spectral index dis¬ 
tribution and significantly under-predicts the number of sources in 
the faintest flux density bin. It is likely that the observed faint, flat- 
spectrum sources are a result of the cores of Fanaroff and Riley type 
I (FRI; Fanaroff & Riley 1974) sources becoming dominant at high 
frequencies. These results highlight that the faint, high-frequency 
source population is poorly understood and therefore the import¬ 
ance of further study of this population. 

In this paper, we select all sources from the IOC sample stud¬ 
ied in Paper I which which have deep 1.4-GHz observations avail¬ 
able (from the Guglielmino et al. 2012, Biggs & Ivison 2006 or 
Owen & Morrison 2008 surveys, see Section 2.1 for full details) 
and have flux densities above the IOC completeness limit in the 
given region. This complete sample of 96 sources is matched to op¬ 
tical and infrared data available in the field. These data enable us to 
distinguish between dilferent source types and estimate photomet¬ 
ric redshifts for the objects. Crucially, by using radio-to-optical ra¬ 
tios we are able to separate radio-loud AGN, which we know dom¬ 
inate the extragalactic radio source population at higher flux dens¬ 
ities (5 15 GHz <: 10 mJy) from radio-quiet AGN and starforming 
galaxies, which are predicted to begin to contribute to the source 
population at lower flux density levels (e.g. Wilman et al. 2008). 
The range of multi-wavelength data available in the Lockman Hole 
is described in Section 2 and the methods used to match these cata¬ 
logues to the IOC catalogue are described in Section 3. Photomet¬ 
ric redshifts are calculated in Section 4 and combined with avail¬ 
able spectroscopic redshifts to produce a final redshift catalogue. 
In Section 5 the radio-to-optical ratio is estimated for all sources in 
this sample. The properties of the sample are discussed in light of 
the redshift information in Section 6. The IOC sample is compared 
to the S^ simulation in Section 7 and to other observational studies 
in Section 8. 


2 DATA USED 
2.1 Radio data 

The work in Paper I is based on a sample of sources selected from 
the IOC survey at 15.7 GHz. The IOC survey was observed with 
the Arcminute Microkelvin Imager (AMI; Zwart et al. 2008) which 
has a resolution of 30 arcsec. Full details of the IOC survey can 
be found in AMI Consortium: Davies et al. (2011) and AMI Con¬ 
sortium: Franzen et al. (2011). This work uses the Lockman Hole 
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Figure 1. Positions of all 296 sources in the IOC Lockman Hole sample. 
The positions of the 96 sources in the sample used in this paper are shown 
as small circles (red in the online version), and the remaining 200 sources 
which are not in this sample are shown as crosses (blue in the online ver¬ 
sion). The sources which fall helow the IOC completeness limit, and are 
therefore not included in the sample used in this paper, as shown as pale 
blue crosses. The IOC complete areas are shown by the rectangles (the large 
rectangles indicate the regions complete to 1 mJy and small rectangles con¬ 
tained within these show the regions complete to 0.5 mJy). The WSRT sur¬ 
vey area is shown hy the large (green) circle and the BI2006 and OM2008 
survey areas are shown hy the (purple) small circle and small square. 


field of the survey, which consists of 4.64 deg^ complete to 1 mJy 
and 1.73 deg^ complete to 0.5 mJy. In Paper I we investigated the 
radio properties of this sample by matching the catalogue to sev¬ 
eral lower-frequency catalogues available in the field; a deep Giant 
Meterwave Radio Telescope (GMRT) survey at 610 MHz (Garn et 
al. 2008, 2010), a Westerbork Synthesis Radio Telescope (WSRT) 
survey at 1.4 GHz (Guglielmino et al. 2012), two deep Very Large 
Array (VLA) surveys at 1.4 GHz by Biggs & Ivison 2006 (BI2006) 
and Owen & Morrison 2008 (OM2008), the National Radio Astro¬ 
nomy Observatory (NRAO) VLA Sky Survey (NVSS; Condon et 
al. 1998) and Faint Images of the Radio Sky at Twenty centimetres 
(FIRST; White et al. 1997). These radio data are summarised in 
Table 1. Lower-frequency counterparts were found for 266 out of 
the 296 sources in the sample, allowing radio spectral indices to be 
calculated for these 266 sources and upper limits to be placed on 
the spectral indices of the remaining 30 sources. 

As many of these lower frequency catalogues have higher res¬ 
olutions than the IOC survey, they also provide more accurate po¬ 
sitions which are vital when searching for multi-wavelength coun¬ 
terparts for the IOC. It is therefore useful to define a complete sub¬ 
sample of sources which have 1.4-GHz data available. For this pur¬ 
pose, we used the BI2006, OM2008 and WSRT surveys as they are 
the deepest in the field. OM2008 and BI2006 have rms noises of 
2.7 and 6.0 pjy/beam respectively and all the IOC sources in these 
fields are detected at 1.4 GHz. The majority of the Guglielmino 
et al. WSRT map has an rms noise of < 15 pJy/beam and parts 
have an rms noise of 11 pJy/beam. This means that it should be 
possible to detect the faintest sources in the complete IOC sample 
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Table 1. Radio catalogues in the Lockman Hole used in Paper I. 


Catalogue 

Reference(s) 

Epoch of observation 

Frequency 

Beam size 

rms noise 




/GHz 

/arcsec 

/mJy 

IOC-shallow j 

1 AMI Consortium: Franzen et al. (2011) 

I AMI Consortium: Davies et al. (2011) 

Aug 2008 - June 2010 

15.7 

30 

0.1 

IOC-deep j 

1 AMI Consortium: Franzen et al. (2011) 

1 AMI Consortium: Davies et al. (2011) 

Aug 2008 - June 2010 

15.7 

30 

0.05 

GMRT ] 

I Garn et al. (2008) 

1 Garn et al. (2010) 

Jul 2004 - Oct 2006 

0.610 

6x5 

0.06 

WSRT 

Guglielmino et al. (2012) 

Dec 2006 - Jun 2007 

1.4 

11x9 

0.011 

OM2008 

Owen & Morrison (2008) 

Dec 2001-Jan 2004 

1.4 

1.6 

0.0027 

OMK2009 

Owen et al. (2009) 

Feb 2006 - Jan 2007 

0.324 

6 

0.07 

BI2006 

Biggs & Ivison (2006) 

Jan 2001 - Mar 2002 

1.4 

1.3 

0.0046 

FIRST 

White et al. (1997) 

1997-2002 

1.4 

5 

0.15 

NVSS 

Condon et al. (1998) 

1997 

1.4 

45 

0.45 
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Figure 2. Spectral index and flux density distributions for sources in the full sample of 296 sources used in Paper I and sample of 96 sources discussed in this 
paper. Left: spectral index, right: 15.7-GHz flux density. Note that the spectral index distribution for the full sample includes 30 sources with upper limits on 
their spectral indices. 


(with 5 15 . 7 GHz > 0.5 mJy) in the WSRT map provided they have 
spectral indices > -1 (where S oc v““), assuming a Scr WSRT 
detection. Sources with spectra which rise as steeply as this are 
very rare so the vast majority of the IOC sources in the WSRT 
map are detectable - in fact, all but one of the sources in the com¬ 
plete IOC catalogue are detected in the WSRT map. We therefore 
define the sample as all IOC sources in the complete catalogue in 
the OM2008, BI2006 or WSRT deep survey areas. This sample 
contains 96 sources and accurate positions and spectral index in¬ 
formation are therefore available for all but one of the sources. The 
positions of the sources in the sample are shown in Fig. 1. This 
sample of 96 sources with deep 1.4 GHz data available is the sub¬ 
ject of this paper. 

The spectral indices and flux density distributions of the full 
sample of 296 sources studied in Paper I and the sub-sample of 96 
sources studied in this paper are plotted in Fig. 2. The spectral index 
distributions of the two samples are relatively similar, although all 
nine of the very steeply rising sources, with < -0.8, in the 
full sample are not included in the sample studied in this paper. 
Seven of these nine sources are below the IOC completeness limits 
(0.5 mJy in the deep regions and 1 mJy in the shallow regions) and 
are not detected at 1.4 GHz so the spectral indices are upper limits 


calculated from the 3cr noise in the WSRT map. The flux density 
distributions of the full sample and sub-sample used in this paper 
are also very similar, although none of the faintest sources in the 
full sample appear in the sample studied here, as this sample only 
contains sources above the IOC completeness limits (0.5 and 1 mJy 
in the deep and shallow regions respectively). 

2.2 The Lockman Hole SERVS Data Fusion 

The Spitzer Extragalactic Representative Volume Survey (SERVS; 
Mauduit et al. 2012) is a warm Spitzer survey which imaged 
~ 18 deg- using the 3.6 and 4.5 pm IRAC bands down to an rms 
noise of 0.4 pjy. The SERVS Data Fusion is a multi-wavelength 
infrared-selected catalogue compiled by Vaccari et al. (in pre¬ 
paration, www.mattiavaccari.net/df) containing most of the 
multi-wavelength public photometry and spectroscopy for SERVS 
sources within the Lockman Hole and other SERVS fields. For 
inclusion in the SERVS Data Fusion, a source must be detec¬ 
ted at either 3.6 or 4.5 pm in the SERVS images, and ancillary 
datasets are matched against the SERVS position using a search 
radius of 1 arcsec. In the Lockman Hole, the SERVS Data Fu¬ 
sion includes optical photometry by Gonzalez-Solares et al. (201 1) 
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(GSll), near-infrared photometry from the United Kingdom In¬ 
frared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS; 
see Lawrence et al. 2007) and mid-infrared and far-infrared pho¬ 
tometry from the Spitzer Wide-Area Infrared Extragalactic survey 
(SWIRE; see Lonsdale et al. 2003). 

The GS11 deep optical data were taken with the Wide Eield 
Camera (WFC) at the Isaac Newton Telescope (INT) and the 
Mosaic-1 camera on the May all 4-meter Telescope at the Kitt-Peak 
National Observatory (KPNO) in g, r, i and z-bands. The average 
magnitude limits in the g, r, i and z-bands are 24.5, 24.0, 23.3 and 
22.0 (AB, 5(T for a point-like object measured in a 2-arcsec aper¬ 
ture). Information about the optical morphology of the objects is in¬ 
cluded in the full published GS 11 catalogue but is not in the SERVS 
Data Fusion catalogue. We therefore matched the full GSll cata¬ 
logue to the Data Fusion catalogue to include this information. A 
match radius of 1 arcsec was used and all of the objects with op¬ 
tical information in the Data Fusion catalogue have matches to the 
GSll catalogue. 

UKIDSS (Lawrence et al. 2007) used the UKIRT Wide Field 
Camera (WFCAM) to map 7500 deg^ in five different surveys in 
J, H and K bands. The Lockman Hole is part of the Deep Ex¬ 
tragalactic Survey (DXS), which has a limiting K magnitude of 
21 (vega). The Data Fusion catalogue contains data from UKIDSS 
Data Release 9, which includes only J and K bands in the Lockman 
Hole. 

SWIRE (Lonsdale et al. 2003) is a wide-field high galactic 
latitude survey covering nearly 50 deg^ in six different fields, one 
of which is the Lockman Hole. These fields have been surveyed 
by the Spitzer Space Telescope using both the Infrared Array Cam¬ 
era (IRAC) and the Multi-Band Imaging Photometer (MIPS) far- 
infrared camera. The Data Fusion catalogue contains data from 
IRAC, which made observations at 3.6, 4.5, 5.8 and 8.0 pm with 
5cr sensitivities of 3.7, 5.4, 48 and 37.8 pJy respectively (the longer 
wavelength MIPS photometry is not included here). 

These catalogues are summarised in Table 2. There are a max¬ 
imum of ten photometric bands available for each source. The data 
fusion catalogue also contains the spectroscopic redshifts from a 
number of different catalogues available in the field. The spectro¬ 
scopic redshifts used in this paper and their references are listed in 
table A2. 


2.3 Fotopoulou et al. photometric redshift catalogue (F12) 

Fotopoulou et al. (2012) (F12) produced a deep photometric red- 
shift catalogue covering 0.5 deg^ contained within the southern 
IOC Lockman Hole field using the Le Phare photometric redshift 
code (Amouts et al. 1999; Ilbert et al. 2006). This catalogue con¬ 
tains 187,611 objects and has up to 21 bands available, ranging 
from far-ultraviolet to mid-infrared. The far-ultraviolet (FUV) and 
near-ultraviolet (NUV) observations used in F12 were made by the 
Galaxy Evolution Explorer (GALEX), with limiting magnitudes of 
24.5 in both bands. At optical wavelengths, they use data from 
the Large Binocular Telescope (LBT), Subaru, and the Sloan Di¬ 
gital Sky Survey (SDSS). The LBT data consists of five bands, 
U, B, V, Y and z', and covers about 0.25 deg^, the Subaru data 
contain E and z'-band observations and SDSS contains data in 
u', g', r', i' and z'-bands. F12 also use the UKIDSS and Spitzer 
data which are used in this work and described in more detail in 
Section 2.2. X-ray observations are available for 0.2 deg^ of the 
field, and the 388 X-ray-detected sources, presumed to be AGN, 
are treated differently in the fitting process. 


2.4 Revised SWIRE Photometric Redshift Catalogue (RR13) 

Rowan-Robinson et al. (2013) (RR13) produced an updated 
SWIRE photometric redshift catalogue, which is a revised ver¬ 
sion of the redshift catalogue produced by Rowan-Robinson et 
al. (2008). The revised catalogue uses the Data Fusion multi¬ 
wavelength catalogue (Section 2.2), which provides deeper op¬ 
tical data and more photometric bands than the catalogue used in 
Rowan-Robinson et al. (2008). The redshifts are estimated using 
a two-pass template method based on six galaxy and three AGN 
templates in the first pass and 11 galaxy and three AGN templates 
in the second. AGNs are identified by their optical morphology - 
only objects which appear point-like are fitted with AGN templates. 
This reduces the risk of catastrophic outliers occurring when nor¬ 
mal galaxies are erroneously fitted with AGN templates, but does 
mean that some AGN will be misclassified as normal galaxies. The 
main difference between the RR13 catalogue and the previous 2008 
version is the treatment of AGNs, as dust torus emission is now 
included in the quasar templates. The RR13 catalogue only con¬ 
tains redshift values for objects which were included in the ori¬ 
ginal SWIRE photometric redshift catalogue (Rowan-Robinson et 
al. 2008), so does not include redshifts for every object in the Data 
Fusion catalogue. There are also a small number of objects in the 
RR13 catalogue which are not in the Data Fusion catalogue as these 
appear in the original SWIRE redshift catalogue. 

RR13 contains redshift estimates for 1,009,607 sources in all 
eight of the SWIRE fields and covers all but a small section of the 
IOC Lockman Hole field. 


3 MATCHING THE CATALOGUES 
3.1 Morphology of the radio sources 

Due to the high density of sources in the Data Fusion catalogue it 
is necessary to take into account the structure of the radio sources 
when matching the catalogues, as there may be several optical ob¬ 
jects within the radio contours of extended sources. 

In order to determine whether or not a source is extended the 
ratio of total flux to peak flux density (C = 5i„t/5peak) was calcu¬ 
lated for all sources with a match in either the FIRST, GMRT or 
WSRT surveys or the 324 MHz VLA survey by Owen et al. (2009) 
(OMK2009). Flux density values from these lower-frequency cata¬ 
logues are used instead of the IOC catalogue as they have a higher 
resolution than the IOC catalogue (which has a beam size of 30 arc- 
sec). FIRST, GMRT and OMK2009 all have a synthesised beam of 
xS arcsec, while WSRT has a larger beam of = 10 arcsec. (OM2008 
and BI2006 have significantly smaller beam sizes, = 1 arcsec, mean¬ 
ing that values of C from these catalogues cannot be directly com¬ 
pared to the other C values so they are not used in this ana¬ 
lysis). FIRST and GMRT cover the whole field while WSRT and 
OMK2009 only cover part of the field. Those sources with a match 
in either FIRST or GMRT were classified as extended if either 
Cgmrt or CpiRST > 1-2, otherwise they were classified as compact 
for the purpose of matching. Sources without a match in FIRST or 
GMRT but with a match to OMK2009 were classified as extended 
if CoMK 2009 >1.2 and compact otherwise. All sources with a match 
in WSRT and not in the other three catalogues used here were clas¬ 
sified as extended for the purpose of matching as the resolution of 
the WSRT map is not high enough to ensure these sources are not 
extended. 

There are six sources which are not classified as extended or 
compact because they do not have a counterpart in any of the four 
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Table 2. A summary of the multi-wavelength information used in this work. 


Survey 

Reference 

Band 

Flux density limit (5 (t) 

SWIRE 

Lonsdale et al. (2003) 

3.6, 4.5, 5.8, 8.0 [.im 

3.7, 5.4, 48, 37.8 pjy 

SERVS 

Mauduit et al. (2012) 

3.6, 4.5 [.tni 

1.3, 1.5 pJy 

UKIDSS 

Lawrence et al. (2007) 

J,K 

21 (AB) 

GSll 

Gonzalez-Solares et al. (201 1) 

g, r, i, z 

24.5, 24.0, 23.3, 22.0 (AB) 



Figure 3. The separation distribution when the radio sources in the IOC 
catalogue and the simulated catalogue are matched to the Data Fusion cata¬ 
logue, taking all matches to the radio sources within 30 arcsec. Note that 
the full IOC sample is matched here. 

catalogues used when classifying the sources but which do have 
a match in OM2008 or BI2006, and therefore an accurate posi¬ 
tion. These sources all have angular sizes less than 3 arcsec in the 
OM2008 or BI2006 catalogues so were considered to be compact 
for this purpose. 

There are therefore 38 sources which are compact and 57 
sources which have been classified as extended for the purpose of 
matching but which may in fact not be significantly extended on 
these angular scales. These two groups of sources are treated sep¬ 
arately when identifying optical matches. 

3.2 Matching the catalogues 

The different catalogues were matched using the topcat* software 
package. The density of sources in the Data Fusion catalogue 
is high compared to the potential error in the IOC source posi¬ 
tions (x! 6 arcsec) so we therefore use the more accurate posi¬ 
tions from the lower frequency radio catalogues (typical error x: 1 
arcsec). When there are several positions available for a source, 
they are used in the following order of preference: FIRST, GMRT, 
BI2006/OM2008, WSRT, IOC. For sources which are resolved into 
multiple components in FIRST or GMRT, the position from the 
IOC catalogue was used instead as this gives a best estimate of the 
centre of the flux. For the three sources which have two separate 
components listed in the original IOC catalogue (see Paper I for 
details), the average of the two IOC positions is used. 

The match radius needs to be chosen carefully to avoid false 

' see: http://www.starlink.ac.uk/topcat/ 


Table 3. A summary of the matches to the Data Fusion catalogue found for 
the 96 IOC sources studied in this paper. 


Description 

No. of sources 

Extended - probable match 

21 

Extended - possible match 

23 

Extended - confused 

8 

Extended - no match 

5 

Compact - match within 2 arcsec 

36 

Compact - no match within 2 arcsec 

2 

IOC position only - no matching attempted 

1 


matches while still maximising the number of real matches. The 
IOC sources were shifted by 0.2 degrees in declination to produce 
a simulated sample of randomly positioned sources. Both this sim¬ 
ulated sample and the real sample were matched to the Data Fusion 
catalogue, and all Data Fusion objects within 30 arcsec of each 
source were noted. The separation between the matches is shown 
in Fig. 3; it is clear that beyond 2 arcsec the number of real and ran¬ 
dom matches becomes comparable. (Note that the full IOC sample 
is used in this plot to provide better statistics.) 

For those sources classified as compact, the nearest match 
within 2 arcsec was accepted. If there was no match within 2 arcsec 
then the source in question was considered to have no optical coun¬ 
terpart. In total, 36 of the compact sources have a match within 2 
arcsec and two do not. 

For sources classified as extended the positions of all the op¬ 
tical sources were plotted on top of the radio contours (GMRT maps 
were used for those sources with a GMRT match, WSRT maps were 
used for the remaining sources). These images were then examined 
and the sources were also assigned one of the following flags: 

(i) probable match - only source within the 3(T radio contours 
(21 sources); 

(ii) possible match - looks likely but there are other sources 
within the 3cr radio contours (23 sources); 

(iii) confused - several sources within the 3cr radio contours so 
cannot identify the correct match (8 sources); 

(iv) no match - no sources within the 3cr radio contours (5 
sources). 

Examples of sources assigned each of these flags are given in Fig. 4. 
Table 3 contains a summary of the number of matches to the Data 
Fusion catalogue and the flags used in the catalogue. 

In summary, we have identified possible counterparts for 80 
out of the 96 sources (83 percent). A table listing the multi¬ 
wavelength counterparts for each source is included in table A1 in 
the appendix. 24 of these 80 sources have a spectroscopic redshift 
available in the Data Fusion catalogue, these redshifts and their ref¬ 
erences are listed in table A2. Fig. 5 shows the separation between 
the radio position of each source and the Data Fusion object as¬ 
sociated with it. The contour plot of the one source in the sample 
without an accurate position available was examined by eye; there 
are no possible counterparts within the radio contours. This source 
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10 55 32 30 28 26 24 22 

RIGHT ASCENSION (J2000) 

Coni peak flux = 2.8202E-03 JY/BEAM 
Levs = 9.90EE-05 * (-4, -2, 2, 8,16, 32, 64, 128, 

256,512, 1024, 2048) 



Coni peak flux = 1.5408E-02 JY/BEAM 
Levs = 7.277E-05 * (-4, -2, 2, 8,16, 32, 64, 128, 
256, 512, 1024,2048) 



Com peak flux = 3.5248E-02 JY/BEAM 
Levs = 1.305E-04 * (-4, -2, 2,8,16, 32,64, 128, 
256,512, 1024, 2048) 



Com peak flux = 2.2704E-03 JY/BEAM 
Levs = 7.141 E-05* (-4, -2, 2, 8,16, 32, 64, 128, 

256, 512, 1024,2048) 


Figure 4. Examples of extended sources assigned each of the four flags. Top left = 1 (probable match), top right = 2 (possible match), bottom left = 3 
(confused) and bottom right = 4 (no match). The large cross (x) marks the position of the radio source in the IOC catalogue and the contours are taken from 
the GMRT image. The smaller pluses (+) mark the positions of objects in the Data Fusion catalogue. 


is therefore included in the group of sources without a match in 
future discussions (giving a total of eight sources without a match, 
and a further eight confused sources). 

SERVS 3.6-|.im and 4.5-[.im images of the eight sources 
without a match were examined by eye, and in one case (source 
10CJ105040+573308) there was a source visible in the SERVS 
images. This source lies close to a very bright ( 53.6 pm = 1.07 x 
10“* jiJy) source so is not included in the SERVS catalogue. The 
3.6 (im-flux density for this source was estimated from the image, 
and although we do not have enough information to calculate a 
photometric redshift for this source, this value is used in Section 5 
when calculating radio-to-infrared ratios. 

The spectral index and flux density distributions of the sources 
with and without a counterpart are compared in Fig. 6 . The spec¬ 
tral index distributions of the two groups of sources (those with 


and without a match) are broadly similar, except that a higher pro¬ 
portion of the very steep (a[ 4 ^ > 0.8) sources are unmatched (32 
percent of sources with a[^ > 0.8 are unmatched compared to 
11 percent of sources with < 0.8). This is probably because 
these very steep sources tend to be significantly extended (all of the 
sources with a > 1 have angular sizes > 20 arcsec. Section 6.4) and 
are therefore more likely to be classified as confused when match¬ 
ing. The 15.7-GHz flux density distributions of the two groups of 
sources are similar, although all six of the brightest IOC sources 
have counterparts in the multi-wavelength catalogue. 


3.3 Matching to other photometric redshift catalogues 

The sample was matched to the F12 and RR13 photometric redshift 
catalogues. For the 80 sources which have a counterpart in Data Fu- 
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Figure 6. A compaiison of the spectral index and 15.7-GHz flux density distributions of sources which have a match in the Data Fusion catalogue and those 
which do not. 



Figure 5. Separation between the best radio position for each source and the 
position of its counterpart in the Data Fusion catalogue. Sources which are 
considered extended for the purposes of matching are shown by ‘+’ (blue 
in the online version), compact sources are shown by ‘x’ (red in the online 
version). The circle indicates a separation of 2 arcsec (the match radius for 
compact sources). 

sion the catalogues were matched using the position from the Data 
Fusion catalogue and the nearest match within 1.5 arcsec is accep¬ 
ted (the shifting procedure described in Section 3.2 was repeated 
for each of the catalogues to choose these match radii). This gave a 
total of 53 matches to RR13 and 20 to F12. 


3.4 Possible matches for confused sources 

There are eight sources which were classified as confused when 
matching to the Data Fusion catalogue. Although it is not possible 
to identify a single counterpart for these sources, some useful in¬ 
formation about their nature can be gained by looking at all pos¬ 
sible counterparts within the radio contours. Therefore, all objects 


Table 4. A summary of matches to multi-wavelength catalogues for IOC 
sources used in this work. 


Catalogue 

Number of matches to IOC sources 

Data Fusion 

80 

Spec 

24 

GSll 

59 

RR13 

53 

F12 

20 


Notes: 

a) Spectroscopic redshift from the Data Fusion catalogue. 


within one tenth of the peak flux in the GMRT sub-image were se¬ 
lected as possible counterparts for each IOC source. For the one 
source without a GMRT image, all objects within the 3cr contour 
in the WSRT map were selected instead. In total, 30 possible coun¬ 
terparts were identified for the eight sources. These are included in 
later discussion. 


4 PHOTOMETRIC REDSHIFT FITTING 

The two photometric redshift catalogues described in Section 2 do 
not contain redshift values for all the sources in this sample with 
multi-wavelength data available, so we performed our own photo¬ 
metric redshift fitting. The publicly available photometric redshift 
code Le Phare^ (Arnouts et al. 1999; Ilbert et al. 2006) was used to 
compute photometric redshifts for the sources in this sample with 
counterparts in the Data Fusion catalogue. The code takes an input 
library of spectral energy distribution (SED) templates, which are 
assumed to represent the SEDs of the observed sample, and shifts 
them to a range of redshift values. These templates are then fitted 
to the photometric data, and a least-squares minimisation is used to 
select the best-fitting SED template for each source. The redshift 
of the best-fitting template is then adopted as the redshift estimate. 
The photometric data used here was from the Data Eusion catalogue 
(see Section 2.2), which has up to ten photometric bands available 
for each source. Data at 3.6 and 4.5 qm are available from both 
SWIRE and SERVS for some sources, in which case values from 

^ http://www.cfht.hawaii.edu amouts/LEPHARE/lephare.html 


© 2015 RAS, MNRAS in press, 1-?? 


































































8 I. H. Whittam et al. 


SERVS were used for the photometric fitting, as these observations 
have a better signal to noise. 

Each source was fitted to two ditferent template libraries, the 
first containing galaxy templates and the second containing AGN 
templates. These two libraries, and the extinction laws applied, 
are the same as those used in El2. The galaxy templates used 
are the library produced by Ilbert et al. (2009). These include 
nine templates generated by Polletta et al. (2007) - three ellipt¬ 
ical galaxy SEDs and six spiral galaxy SEDs (SO, Sa, Sb, Sc, Sd, 
Sdm) - and 12 starburst galaxy SEDs generated using the Bruzual 
& Chariot (2003) models (with starburst ages ranging from 3 to 
0.03 Gyr). This gives a total of 21 SED templates. Ilbert et al. 
(2009) linearly interpolated between some of the Polletta et al. 
templates to refine the sampling in colour-redshift space, result¬ 
ing in a total of 31 templates. For templates Sb to SB3 (template 
IDs 11 to 23) extinction is applied according to the Prevot et al. 
( 1984) Small Magellanic Cloud (SMC) law, while for SB4 to SB 11 
(template IDs 24 to 31) the Calzetti et al. (2000) laws are ap¬ 
plied. No additional extinction is applied for templates earlier than 
Sb. The intrinsic galactic absorption is calculated with values of 
E{B-V) = 0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.40 and 0.50. 
Emission lines were added to the templates using the option in the 
Le Phare code as this has been shown to give better results, even in 
the case of broadband photometry (Ilbert et al. 2009). 

The AGN template library was taken from Salvato et al. 
(2009). This library of 30 templates contains galaxy and AGN tem¬ 
plates, as well as a number of hybrid templates. These hybrid tem¬ 
plates contain contributions from both galaxy and AGN templates 
in proportions ranging from 10%: 90% to 90%: 10% in steps of 10% 
(see Salvato et al. 2009 for full details). For this template library 
extinction is allowed to vary from E(B - V) = 0 to 0.5 in steps of 
0.05. 

For both libraries the templates are calculated at redshifts 0-6 
in steps of Az = 0.01 and in steps of Az = 0.02 for redshifts 6-7. 
The absolute magnitude was restricted to the range -28 < M < -8, 
where M is the absolute magnitude in the Al-band. 

4.1 Selecting a final redshift value 

Each of the 80 matched objects in the sample has been fitted to two 
SED template libraries (the galaxy template library and the AGN 
template library), so there are two possible redshift values produced 
for each source. Throughout this section the redshift value resulting 
from the best-fitting template from the galaxy library is referred to 
as Zgal and the redshift value resulting from the AGN template 
library is referred to as zagn- 

To decide which value is the most appropriate one to use for 
each source the SEDs which provided the best fit to the photometric 
data from each of the two libraries were plotted for each source 
and examined by eye. Both the galaxy and AGN template fits were 
then qualitatively assigned one of the following template flags to 
characterise the fit: 

1 = good fit; 

2 = possible fit; 

3 = poor fit. 

One of the possible two redshift values was then selected for each 
source, zgal was selected if the fit to the galaxy template is bet¬ 
ter than the fit to the AGN template (i.e. galaxy template flag < 
AGN template flag). If the fit to the galaxy and AGN templates were 
judged to be equally good (i.e. galaxy template flag = AGN tem¬ 
plate flag), then zgal was selected if the object’s optical morpho¬ 
logy indicated that it was extended in the optical observations, and 


Zgal was also selected if there was no optical morphology inform¬ 
ation available, zagn was selected if the fit to the AGN template 
was better than the fit to the galaxy template (i.e. AGN template fit¬ 
ting flag < galaxy template fitting flag), or if the two fits appeared 
equally good and the object appeared point-like in the optical. Any 
source with fewer than three photometric bands available was not 
assigned a redshift value. 

In the final photometric redshift catalogue 59 sources have 
redshift values from the galaxy templates library and 11 have red¬ 
shift values from the AGN templates library. This gives a total of 70 
sources with redshift values from the Le Phare photometric redshift 
fitting process (the remaining ten sources with matches to the Data 
Fusion catalogue had fewer than four photometric band available 
so no redshift estimate is included for these sources). 

4.2 Comparison with other redshift catalogues 

The redshift values obtained from the Le Phare photometric red¬ 
shift fitting were compared to the F12 and RR13 photometric red¬ 
shift catalogues. First the two published catalogues were compared 
to each other. The two catalogues were matched using a match 
radius of 1.5 arcsec, giving 7895 matches; the redshift values of 
these matched sources from the two catalogues are compared in 
Fig. 7. This plot shows significant scatter, demonstrating that there 
are some significant disagreements between the two catalogues and 
highlighting the difficulties in achieving reliable redshifts using 
photometric methods for these sources. This is demonstrated fur¬ 
ther in Fig. 8 which shows the percentage difference between the 
redshift values from the two catalogues. The percentage of ‘cata¬ 
strophic’ outliers, defined as when Z 2 - Zi/(1 + z) > 0.3, where 
z = (zi -t- Z 2 )/ 2 , is 20 per cent. Redshift values derived in this 
work are compared to the F12 and RR13 catalogues in Fig. 9. For 
the majority of the sources there is a good agreement between the 
Le Phare values calculated here and the values from RR13 and 
F12. However, for several sources there are significant differences 
between the Le Phare values and the values from the other two 
catalogues, with the percentage of catastrophic outliers being 24 
per cent when compared to the El2 catalogue and 30 per cent when 
compared to the RR13 catalogue. Many of the sources do not agree 
within the error bars, although note that these error bars simply 
quantify the goodness of the fit of the photometric data to the 
chosen template, which does not necessarily mean that the chosen 
template is the correct one. The right panel of Fig. 9 shows the Le 
Phare photometric redshifts compared to the spectroscopic values 
for the 24 sources with spectroscopic values available. There is a 
reasonable agreement for the majority of the sources, with a cata¬ 
strophic outlier percentage of 29 per cent. 

4.3 Compiling the final redshift catalogue 

There are spectroscopic redshifts available for 24 sources from the 
Data Fusion catalogue, which are listed in table A2 in the appendix. 
For the remaining sources a photometric redshift is selescted in the 
following order of preference: 1) photometric redshift from RR13 
(38 sources), 2) photometric redshift from F12 (6 sources) and 3) 
photometric redshift from Le Phare fitting (9 sources). RR13 val¬ 
ues were given preference over El2 values because they are avail¬ 
able for a greater number of sources as the RR13 catalogue covers 
the whole IOC survey area, giving greater consistency. A full cata¬ 
logue of the redshifts is given in table Al in the appendix. 

The redshift distribution of all 77 sources with a redshift es¬ 
timate (or value) is shown in Fig. 10. The median redshift is 0.91 
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Figure 7. Comparison of the redshift values from the RR13 and F12 catalogues. The dashed line indicates where the redshift values from the two catalogues 
are equal. The same values are plotted in the two panels - in the left-hand panel log(l + z) is plotted to make it easier to see the high density of points, and the 
right-hand panel shows z on a linear scale for easier comparison with Fig. 9. 




^RR13 ^F12 ^spBC 


Figure 9. Comparison of the redshift values from our work using Le Phare with those from the RR13 (left panel) and F12 (middle panel) catalogues and the 
spectroscopic redshift values (right panel). The dashed line indicates where the redshift values from the two catalogues are equal. Icr error bars are plotted for 
sources in the Le Phare catalogue with a final value from the galaxy library and values from the F12 catalogue. Sources with only three photometric bands 
available for the Le Phare fitting are circled (in red in the online version); these values are not included in the final catalogue. 



Figure 8. Percentage difference between the redshift values from the RR13 
and F12 catalogues (i.e. 100 X (zrro - zfi2)/(1 + z) where z = (zrrb + 
4F12)/2). 


with an interquartile range of 0.84. The redshift distribution of the 
sample is discussed in Section 6. 

4.4 Confused sources 

In Section 3.4, a total of 30 possible counterparts were identified 
for the eight ‘confused’ sources. The Le Phare photometric code 
was run on these possible counterparts in exactly the same way 
as detailed in Section 4. The resulting redshift values for all the 
possible counterparts to each ‘confused’ IOC source are shown in 
Fig. 11. 


5 RADIO TO OPTICAL RATIO 

The ratio of the radio and the optical flux densities of a source 
provides useful information about its nature, since for the same op¬ 
tical magnitude radio galaxies have much higher radio flux densit¬ 
ies than either radio-quiet AGN or starforming galaxies. The radio- 
to-optical ratio R was defined following Condon (1980), using: 

= Sl.4GHzX lO"""'"-*"-’’ (1) 
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Figure 10. Redshift distribution for all sources with a redshift estimate. 
Sources with spectroscopic redshift values are shown separately. 
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Figure 11. Redshift values for the possible counterparts for each IOC source 
classified as confused. 


where S 1.4 ghz is the flux density at 1.4 GHz in mJy and m is the 
optical magnitude in the i-band. Sources with radio-to-optical ra¬ 
tios R > 1000 are considered to be radio loud (Machalski & Con¬ 
don 1999), while those with smaller values of R are classified as 
radio quiet (and are therefore either radio-quiet AGN or starform¬ 
ing sources), f-band magnitudes are used here, as they are the best 
match to the /-band magnitudes used by Machalski & Condon. 

The radio-to-optical ratio, R, was calculated using equation 1 
for the sources which have an f-band magnitude available. Of the 
80 sources with a match in the Data Fusion catalogue, 44 have an 
i-band magnitude, whilst 36 sources are not detected in the i-band. 
For the sources without an i-band magnitude, including the eight 
sources which are unmatched in the Data Fusion catalogue, lower 
limits on R are calculated using a lower limit on the i-band mag¬ 
nitude of 23.3 (the limiting magnitude of the i-band observations). 
For the eight sources classified as confused when matching, R was 


Figure 12. The distribution of radio-to-optical ratio, R. Sources with an i- 
band detection are in black. Lower limits are included for those sources 
without an i-band detection (grey) and for those classified as confused 
(white); these sources could move to the right on this diagram. Sources 
with R > 1000, indicated by the vertical dashed line, are considered radio 
loud. 


calculated using the brightest of the possible counterparts identified 
in Section 3.4, which serves as a lower limit on R. 

The distribution of the radio-to-optical ratio of all sources is 
shown in Fig. 12, including lower limits for unmatched and con¬ 
fused sources. Based on the study out to z ~ 0.2, Machalski & Con¬ 
don (1999) found that 98 percent of normal and starburst galaxies 
had R < 1000, so this is used as the cut-off point between radio- 
loud and radio-quiet objects. Fig. 12 shows that only three sources 
have R < 1000, so using this criterion the other 93 are radio loud. 

There are far more detections in the SERVS mid-infrared 
bands than in the optical i-band traditionally used to calculate 
radio-to-optical ratios. It is therefore useful to calculate a ‘radio- 
to-optical’ ratio based on 3.6-|.im flux densities (referred to as the 
radio-to-infrared ratio, R^ f, ^n,, from now on) from SERVS as there 
are far fewer lower limits. There is the additional advantage that 
the 3.6-pm band remains longwards of the 4000 A break out to red- 
shifts of z ~ 8 . By contrast, the i-band samples below the break at 
z 2: 1 ; variations in starformation rate and absorption by dust have 
strong effects in that part of the spectrum and can introduce uncer¬ 
tainties in interpreting the R values based on i-band magnitudes at 
higher redshifts (Seymour et al. 2008). Calculating 3.6-pm radio- 
to-infrared ratios therefore also enables us to check that this is not 
having a major effect on our results. 

We therefore define a radio-to-infrared ratio based on 3.6-|im 
flux densities (/?3 5 ^m) as follows: 


R 


3.6 [.im 


S 1.4 GHz 
‘^3.6 ;im 


( 2 ) 


80 sources have a 3.6-[rm flux density available, including one 
source where the flux density is not included in the SERVS cata¬ 
logue and is estimated from the map (see Section 3.2). Lower lim¬ 
its on /?3 6 nm were calculated for the eight sources without a 3.6-|xm 
band detection using the SERVS 3.6-pm 5 cr detection limit as an 
upper limit on the 3.6-pm flux density of the sources. Lower limits 
on Ri f, were calculated for the eight confused sources using the 
highest 3.6-pm flux density of any of the possible counterparts for 
each confused source. 
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Figure 13. Radio-to-optical and infrared ratios for !-band and 3.6 pm re¬ 
spectively. All sources with values of both R and R 3.6 |im available ai'e shown 
as (blue) those which are at redshifts of z < 1 are circled (in red in 
the online version). The line of best fit to all sources is shown as a (blue) 
solid line, the line of best fit to the sources with z < 1 is shown as a (red) 
dashed line. The vertical dotted line is at R = 1000, the value used to dis¬ 
tinguish between radio-loud and radio-quiet sources. The horizontal dotted 
line shows the R 3 5 dividing value derived from the best fit line to sources 
with z < 1 and the horizontal dot-dashed line shows the R 3 j cut-otf de¬ 
rived from the best fit line to all sources with both R values available. 



Figure 14. The distribution of radio-to-infrared ratio, R 3 sum- Sources with 
a 3.6 pm-band detection are in black. Lower limits are included for those 
sources without an 3.6 pm-band detection (grey) and for those classified as 
confused (white); these sources could move to the right on this diagram. 
Sources with R 3 6nm > 3.1, indicated by the vertical dashed line, are con¬ 
sidered radio loud. 


Figure 15. Radio-to-optical and radio-to-infrared ratios for all sources. 
Sources which only an infra-red R 3 5 value have a lower limit on R and 
could therefore move to the right. Sources which only have an optical R 
value could move up. Sources with a lower limit on both values of R could 
move up or to the right. The dashed lines show the values of R used to 
distinguish between radio-quiet and radio-loud sources. 


sampling below the 4000 A break in the t-band, were not having a 
significant effect on this choice of cut-off the analysis was repeated 
using sources with redshifts of z < 1 only. The redshifts estimated 
in Section 4 were used to perform the redshift cut, which resulted in 
20 sources being excluded from the analysis. Both fits and the cut¬ 
offs derived from them are shown in Fig. 13. Excluding the sources 
with z > 1 resulted in a lower cut-off of R 3 5 = 3.1. This value is 

used to distinguish between radio-quiet and radio-loud sources in 
the following analysis. A histogram of the R 3.6 ixm values is shown 
in Fig. 14. 

Fig. 15 shows the radio-to-optical and radio-to-infrared ratios 
for all 96 sources in the sample, including lower limits where ne¬ 
cessary. It is clear from this figure that only one source is signific¬ 
antly radio quiet, while several other sources lie close to the radio- 
loud/radio-quiet boundary. There are six sources which could be 
classified as radio quiet using one of the two ratios; three sources 
are classified as radio quiet using the R 3 j criterion only and three 
are radio quiet using both criteria (one of these classifications is 
based on lower limits on both ratios so the source could actually 
be radio loud). Therefore at least 90 out of the 96 sources (94 per 
cent) in this sample are radio loud; in fact, given how close to the 
boundary five of the remaining six sources are, it is likely that they 
are also dominated by AGN activity. This backs up the conclusions 
of Whittam et al. (2014), in which we used VLBI data to show that 
at least 65 per cent of the IOC sources are associated with an AGN. 
Thus at ~ 1 mjy the faint, high-frequency sky is still dominated by 
radio-loud AGN rather than by starforming galaxies or radio-quiet 
AGN. 


The sources with both R and R 3 5 values available (43 
sources) were used to define the value of R 3 ^ used as the cut-off 
between radio-quiet and radio-loud sources. A linear best fit line 
was fitted to the data, and this was used to estimate an cut¬ 

off equivalent lo R = 1000, as shown in Fig. 13. This gave a cut-off 
value of R 3 6 = 4.1. In order to check that the sources with red¬ 

shifts greater than z ~ 1 , where there could be uncertainties due to 


5.1 Far-infrared - radio correlation 

Radio-loudness can also be defined in terms of the far-infrared - ra¬ 
dio correlation (e.g. Ibar et al. 2008; Luchsinger et al. 2015) which 
provides a useful comparison as it gives a view of radio loudness 
not affected by the obscuration of the AGN in the optical. The far- 
infrared - radio correlation is often quantified using the far-infrared 
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Figure 16. Far-infrared - radio correlation parameter qiR. Upper limits are 
shown in white and could move to the left. 


three sources classified as radio quiet based on their radio-to-optical 
ratios and the six sources classified as radio quiet based on their 
radio-to-infrared ratios are circled in black on the left and right pan¬ 
els of Fig. 18 respectively. As for the ratios determined at 1.4 GHz, 
only one of these radio-quiet sources has 15.7-GHz radio-to-optical 
and radio-to-infrared ratios significantly lower than the rest of the 
sample; the other radio-quiet sources all have 15.7-GHz ratios sim¬ 
ilar to the rest of the sample. This significantly radio-quiet source 
(IOC J105028-I-574522) has the lowest redshift in the sample, with 
z = 0.072 (spectroscopic redshift) and has a rising spectra. 

This analysis shows that all of the radio-quiet sources have 
rising spectra. They are unlikely to be starforming galaxies, as star¬ 
forming galaxies only have rising spectra at 15 GHz at z > 3 
(Murphy et al. 2014), and although our data is very deep it is 
not deep enough to detect starforming galaxies at z > 3 unless 
they have a prodigious amount of star formation. These sources are 
therefore likely to be dominated by AGN emission. 


- radio ratio qiR (e.g. Appleton et al. 2004; Garn & Alexander 
2009), where 

qiR = loglo(S24 nm/Sl.4 GHz) (3) 

Starforming sources and radio-quiet AGN are expected to have pos¬ 
itive qiR values, for example, Marleau et al. (2007) found that the 
typical value of q/R for starforming galaxies was 0.83±0.31. Radio- 
loud objects tend to have much lower values of qiR, with typical 
values -0.6 to -1.2 (Prandoni 2010b). 

23 sources have a 24 pm-flux density available from the 
SWIRE catalogue and these were used to calculate qiR values, 
which are shown in Fig. 16. For the undetected sources, the 5cr 
detection limit of 450 pjy was used to calculate an upper limit on 
qiR. Fig. 16 shows that the majority of the sources have qiR val¬ 
ues consistent with those expected for radio-loud AGN. Only eight 
sources have qiR > 0, four of which are upper limits so could have 
much lower values, and only one source has qiR > 0.5 (this is the 
same source that was the the only source to be identified as signific¬ 
antly radio quiet in Fig. 15). This provides further evidence that this 
sample is entirely dominated by radio-loud sources, and shows that 
this result is not affected by obscuration of the AGN in the optical. 


6 RADIO SOURCE PROPERTIES 

In this section, the radio properties of the sources are considered 
in light of the redshift values derived in Section 4. Redshift val¬ 
ues are available for 78 sources, nearly a third of which are spec¬ 
troscopic reshifts. The large errors on some of the redshift estim¬ 
ates, along with the significant discrepancies between catalogues in 
some cases, mean that photometric methods cannot be used to pro¬ 
duce a reliable redshift for any individual source. They can, how¬ 
ever, be used to provide valuable information about the properties 
of the population as a whole. 

6.1 Summary of the radio properties of the sample 

The radio properties of the sources in this sample (along with those 
of the full IOC Lockman Hole sample) were presented in Paper I. 
Flux densities are available at a range of frequencies including 
1.4 and 15.7 GHz, and radio spectral indices are calculated for all 
sources but one (an upper limit is available for the one remain¬ 
ing source). The sources were split into flat-spectrum sources, with 
a < 0.5, and steep-spectrum sources with a > 0.5. 


5.2 Correlations between R and spectral index 


Fig. 17 shows the radio-to-optical and radio-to-infrared ratios as 
a function of radio spectral index. There appears to be a positive 
correlation between both ratios and spectral index; this is, however, 
a selection effect as the sample is selected at 15.7 GHz but 1.4-GHz 
flux densities are used to calculate the ratios. To test for any real 
correlation between the ratios and spectral index radio-to-optical 
and radio-to-infrared ratios were calculated using 15.7-GHz flux 
density, as follows. 


Ri5 - S 1 


. X 10' 


|0.4(«7-12.5) 


R 


3.6 [im 15 ~ 


S 15.7 GHz 
5 3.6 |im 


(4) 


These ratios are shown as a function of spectral index in Fig. 18. 
The positive correlation is no longer seen, confirming that this ef¬ 
fect was due to the sources being selected at 15.7 GHz. 

Fig. 17 shows that all six sources classified as radio quiet us¬ 
ing either or both ratios have rising spectra (a < 0). This rising 
spectral shape means these sources have lower 1.4-GHz flux dens¬ 
ities than the majority of the IOC sample, which could explain why 
they have small radio-to-optical and radio-to-infrared ratios. The 


6.2 Redshift distribution 

There are redshift values or estimates for 77 out of the 96 sources 
in this sample (Section 4). The redshift distributions of steep and 
flat spectrum sources are shown separately in Fig. 19; the redshift 
distributions for the two samples appear fairly similar, although the 
distribution for flat spectrum sources peaks at a slightly higher red¬ 
shift. A KS test was performed on the two samples and the prob¬ 
ability of them being drawn from the same population was 0.16, 
indicating that the two distributions are not significantly different. 


6.3 Luminosity distribution 

Luminosities were calculated for all sources with a redshift value 
or estimate. The luminosities were fe-corrected based on their radio 
spectral index using the following expression: 

U=47rdR(zfSA(l+zr-‘] (5) 

where di is the luminosity distance. 
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Figure 17. Radio-to-optical and radio-to-infrared ratios as a function of spectral index. Left-hand panel shows R and right hand panel shows R 3 g Triangles 
(red in the online version) indicate lower limits on the ratio. The dashed line indicates the cut-off used to distinguish between radio-loud and radio-quiet 
sources (R = 1000 in the left panel and R 3 5 |.im = 3.1 in the right panel). 




Figure 18. Radio-to-optical and radio-to-infrared ratios calculated using 15.7-GHz flux densities as a function of spectral index. Left-hand panel shows R 15 
and right hand panel shows R 3 5 15 - Triangles (red in the online version) indicate lower limits on ratios. The three sources which are classified as radio-quiet 

based on their R value are circled in the left panel and the six sources classified as radio-quiet based on their R 3 5 values are circled in the right-hand panel. 



Redshift 


Figure 19. Redshift distribution for all sources with a redshift value or es¬ 
timate, with steep and flat spectrum sources shown separately. 


The 15.7-GHz luminosity distribution for all these sources is 
shown in Fig. 20, with steep and flat spectrum sources shown sep¬ 
arately. The distributions for the steep and flat-spectrum sources 
are very similar, suggesting there is no difference in luminosity 
between the two populations. The sources display a large range of 
15.7-GFlz luminosities, comparable with those of the powerful FRI 
and FRII radio galaxies (FanarolF & Riley 1974). 


6.4 Linear size distribution 

To calculate linear sizes we need to first estimate angular sizes 
for the sources. Angular size information is available on a range 
of scales from the different radio catalogues available in the field. 
Best estimates of the angular sizes are compiled from these cata¬ 
logues for the 96 sources in this sample. The two deep VLA sur¬ 
veys, BI2006 and OM2008, are the highest resolution surveys, with 
resolutions of 1.3 and 1.6 arcsec respectively, so provide informa¬ 
tion on the smallest angular scales. However, the angular sizes from 
these two catalogues may not be reliable for significantly extended 
sources, as some structure may be resolved out due to a lack of 
short baselines. Additionally, some sources may be resolved into 
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Figure 21. Angular and linear size distributions. The left panel shows angular sizes for all sources in the sample and the right panel shows linear sizes for the 
77 sources with a redshift available. Upper limits are shown in white, and could move to the left. 



log(15.7-GHz luminosity / W Hz^^) 

Figure 20. 15.7-GHz luminosity distribution for all sources with a redshift 
estimate; the sample is divided into fiat and steep spectrum sources. 


multiple components which are then listed separately in the cata¬ 
logue, so the size listed in the catalogue would be for only part 
of the source. For this reason, for any source which was classi¬ 
fied as extended in Paper I the angular size was estimated from the 
Gam et al. (2008) GMRT observations instead. The angular sizes 
for these sources were measured by hand, to avoid any problems 
caused by fitting Gaussians to extended sources with complicated 
structures. For the remaining sources (which were not classified 
as extended) the angular size was taken from the catalogue with 
the highest resolution available. The maximum angular size from 
the relevant catalogue is used, with the catalogues being used in 
the following order of preference: OM2008/BI2006, FIRST, WSRT 
(highest to lowest resolution). A value is flagged as an upper limit 
on the angular size if the size listed is less than the synthesised 
beam size for those observations. 


Two sources required different treatment. The very extended 
source 10CJ105437-I-565922 is not detected in FIRST or GMRT as 
it has diffuse low-brightness structure; its size is measured from the 
NVSS map. The other source, 10CJ104927-I-583830, is only in the 
IOC catalogue, so for this source the IOC beam size (30 arcsec) is 
used as an upper limit on the angular size of the source. 

Linear sizes were then calculated from these angular sizes for 
the 78 sources with redshift information using: 


where D is the linear size of the source, d/, is the luminosity dis¬ 
tance and 6 is the angular size of the source in radians. The angular 
and linear size distributions are shown in Fig. 21; the left panel 
shows the angular sizes and the right panel shows the linear sizes. 
Note that the 18 sources without a redshift value are not included 
in the linear size plot. 

Fig. 22 shows the linear size distributions for flat and steep 
spectum sources. It is evident that the flat spectrum sources are on 
average smaller than the steep spectrum sources. This is expected 
as the extended lobes present in many radio galaxies have steep 
spectra due to optically-thin synchrotron emission, while the cores 
of radio galaxies generally have flat spectra due to self-absorption. 
The spectral index of a source therefore informs us about the relat¬ 
ive contributions of the cores and lobes to the total flux density of 
the source; if a source has a steep spectrum its emission is domin¬ 
ated by the lobes, and it is therefore more likely to appear extended. 
If, however, a source has a flat spectrum it is dominated by emis¬ 
sion from its core, so its lobes may be very weak or not visible at 
all, and the source appears more compact. 

While this trend is tme for the majority of the population, both 
steep and flat spectrum sources cover the full range of linear sizes 
(I ^ D < 1000 kpc), so there are a number of small steep spectrum 
sources and large flat spectrum sources. This is consistent with the 
findings in Paper I and Whittam et al. (2014) which indicate that 
there are small populations of both extended, flat spectrum sources 
and compact, steep spectrum sources in the IOC sample. 

Fig. 23 shows linear size as a function of redshift for all 
sources in the sample used in this paper. This shows that the sources 
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Figure 22. Linear size distribution for sources with redshift values. The top 
panel shows flat spectrum sources and the bottom panel shows sources with 
steep spectra. Upper limits are shown in white, and could move to the left 
on these plots. 



Figure 23. Linear size as a function of redshift. Triangles (red in the online 
version) are upper limits on size, crosses (blue in the online version) are 
values. 


with larger linear sizes tend to be a lower redshifts, and at z > 2 all 
but two sources are unresolved. 


7 COMPAmSON WITH THE SKADS SIMULATED SKY 

In Paper I a sample of sources with 5 ig ghz > 0-5 mjy was selected 
from the catalogue; this sample should be directly comparable 
with the IOC sample. The radio properties of this sample were 
compared to the IOC sample, which showed that the simulation 
fails to accurately reproduce the spectral index distribution of the 
observed sample. The number of flat spectrum sources is massively 



Figure 24. Normalised redshift distribution for sources in and in the IOC 
sample. Objects with no redshift value have been omitted, but included in 
the normalisation. 


underpredicted; there are essentially no sources in the simulated 
sample with a < 0.3, while 40 percent of the IOC sample have 

<03 
“o.6i ^ 

The simulation predicts that the IOC sample should be domin¬ 
ated by FRI sources, making up 71 percent of the population, while 
FRII sources are the second-largest source type (13 per cent). Star¬ 
forming sources make up seven per cent of the simulated sample. 
There is some indication that this may be an overestimate, as we 
find no starforming sources in the IOC sample. However, we are 
sampling a small area so cosmic variance is very high at low red- 
shifts and could therefore be the cause of this discrepancy. 

Fig. 24 shows the normalised redshift distribution for the 
and IOC samples. The IOC sample is normalised by the total 
sample size (96 sources), which includes the 19 sources with 
no redshift value available. The redshift distributions of the two 
samples are similar, although the observed sample displays a 
sharper peak and is missing the extensive low redshift tail in the 
simulated sample. Eleven percent of the sources in the simulated 
sample have z < 0.2, so we would expect to find approximately ten 
sources in the observed sample in this range if the two distributions 
are similar. However, there are only three sources in the observed 
sample with a redshift less than 0.2. The majority of the sources 
with z < 0.2 in the simulation are starforming sources, so the lack 
of low-redshift sources in the IOC sample is consistent with the fact 
that there are very few starforming sources in this sample. 

It is also a possibility that some of the 18 sources which lack 
redshift information and are therefore missing from the observed 
sample have redshifts less than 0.2 and are therefore responsible for 
the discrepancy in the redshift distributions. It is unlikely, however, 
that any of the eight sources without a match are at low redshift, 
as they would have to be very faint in the optical to have z < 0.2 
and not be detected in the optical or infrared observations. Only 
two of the 30 possible counterparts for the confused sources have 
z < 0.2, so these sources cannot account for the missing sources in 
this redshift range. 

The simulation therefore incorrectly predicts that there is a 
population of low-redshift starforming galaxies in the IOC sample. 
This could indicate that the spectra assumed for the starforming 
sources in the simulation is wrong, and they in fact have much 
steeper spectra, or that the extrapolation of the luminosity func- 
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tion for the starforming galaxies is not correct. Studies of the faint 
(5 1.4 GHz <0.1 mjy) source population at lower frequencies by 
Simpson et al. (2012), Lindsay et al. (2014) and Luchsinger et al. 
(2015) have also found fewer starforming galaxies than predicted 
by the simulation. These results support our suggestion that the 
faint end of the luminosity function for starforming galaxies in the 
simulation is not correct, with the number of starforming galaxies 
being overestimated. 

The peak in the redshift distribution also appears to be shif¬ 
ted to slightly lower redshifts in the observed sample, with more 
sources in the bins 0.56 < z < 1.26 (-0.25 < log(z) < 0.1) but 
fewer sources in the bins z > 1.8 (log(z) > 0.25). It is plausible that 
some of the sources without an optical counterpart have z > 1.8, 
and are therefore the cause of this difference. 


8 COMPARISON WITH OTHER STUDIES 

Mignano et al. (2008) investigated the properties of a complete 
sample of 131 radio sources with S > 0.4 mJy observed at 1.4 
and 5 GHz as part of the Australia Telescope ESO Slice Project 
(ATESP) 5 GHz radio survey (Prandoni et al. 2006). This sample 
provides a useful comparison as it has a comparable flux dens¬ 
ity limit to the IOC survey, albeit at a lower frequency. The ESO 
Deep Public Survey provides deep multi-colour (UBVRIJK) im¬ 
ages which cover most of this field, and optical/near-infrared coun¬ 
terparts are found for 66 out of the 85 (78 percent) sources in the 
area covered. Estimates of redshift and optical object type are ob¬ 
tained for 56 of these 66 sources. These results showed that 78 per¬ 
cent of the ATESP 5 GHz sample with optical identifications had 
an active nucleus (i.e. they are either quasars or radio galaxies asso¬ 
ciated with early-type objects), significantly lower that the propor¬ 
tion of radio galaxies found in the IOC sample (94 percent). This is 
confirmed by looking at the radio-to-optical light ratios of the two 
samples; approximately 30 percent of the full ATESP sample (131 
sources) have R < 1000 (classifying them as radio quiet), com¬ 
pared to just six percent of the IOC sources. This suggests selection 
at a higher frequency at this flux density level preferentially selects 
radio galaxies, as the steep-spectrum starforming galaxies drop out 
of the sample. 

Mignano et al. find that those sources in the ATESP sample 
which have flat or inverted radio spectra and are associated with 
objects with early-type spectra are preferentially compact (with lin¬ 
ear sizes < 10 to 30 kpc). They suggest that these sources may be 
ERIs, due to their low radio powers (P 1.4 ghz ~ 10^^“^^^ WHz^') 
and the absence of emission lines in their optical spectra. They do, 
however, note that they would expect ERI sources to have larger 
linear sizes and steeper spectra. As these sources have flat spectra, 
we would expect them to be present in significant numbers in the 
15-GHz-selected IOC sample; these are likely to be the flat, core- 
dominated radio-loud sources which we observe in the IOC sample 
(but which are not present in, for example, the simulation). 

Prandoni et al. (2010a) followed up a sample of early-type 
galaxies selected from the ATESP 5 GHz survey at 4.8, 8.6 and 

19 GHz to further investigate their properties. The main aim was 
to establish whether the AGN population of the sub-mjy sample 
is more closely related to efficiently accreting systems (such as 
radio-quiet quasars), or to systems with low accretion rates (such as 
ERI galaxies), or to low radiative efficiency accretion flows. They 
compare this AGN population to the much brighter (> 500 mJy) 

20 GHz AT20G Bright Source Sample (Massardi et al. 2008) and 
find strong similarities in the radio spectra of the two samples. They 


therefore conclude that the ATESP AGN sources are lower lumin¬ 
osity counterparts of the AT20G ERII radio galaxies, and do not 
find any compelling evidence for a radio-quiet AGN population. 
This is consistent with the properties of the IOC sources, as the 
vast majority of the sources are radio loud. 


9 CONCLUSIONS 

In this paper we have studied the multi-wavelength properties of a 
complete sample of 96 sources selected from the IOC survey. This 
sample was matched to the Data Fusion multi-wavelength cata¬ 
logue, which contains up to ten photometric bands in the optical 
and mid-infrared, and counterparts were identified for 80 out of 96 
IOC sources. Spectroscopic redshifts are available for 24 sources. 

Photometric redshifts were estimated for all sources with suf¬ 
ficient photometric information available using the Le Phare code. 
This produced redshift estimates for 70 of the 80 sources, albeit 
with large errors in some cases. The results are compared to two 
published photometric redshift catalogues (F12 and RR13), and 
are generally in good agreement, although there are some signi¬ 
ficant outliers. These catalogues were then combined to produce a 
final redshift catalogue, which contains redshift estimates for 77 
(24 spectroscopic and 54 photometric) out of the 96 sources in 
the sample. The median redshift of the sample is 0.91 with an in¬ 
terquartile range of 0.84. The large errors on some of the redshift 
estimates, along with the significant discrepancies between cata¬ 
logues in some cases, mean that photometric methods cannot be 
used to produce a reliable redshift for any individual source. They 
can, however, provide information about the properties of a popu¬ 
lation as a whole. 

The radio-to-optical and radio-to-infrared ratios (or lower lim¬ 
its) were calculated for all 96 sources using i-band magnitudes and 
3.6-qm flux densities respectively. Using these ratios six sources 
are classified as radio quiet using at least one of the two values. 
Only one of these sources is significantly radio quiet; the other five 
lie close to the radio-loud/radio-quiet boundary. Therefore at least 
90 out of the 96 sources (94 per cent) in the sample are radio loud, 
indicating that the IOC sample is dominated by radio galaxies. All 
six potentially radio-quiet sources have rising spectra, ruling out 
the possibility that they are starforming galaxies, so their emission 
is dominated by AGN activity. 

These results confirm the conclusions of Paper I that the faint, 
flat-spectrum sources which are found to dominate the IOC sample 
below ~ 1 mJy are the cores of radio galaxies. 

The overall radio properties of the sources in the sample are 
discussed in light of this redshift information; luminosities and lin¬ 
ear sizes are derived for those sources with redshift estimates. There 
is a large range of 15.7-GHz luminosities, with values comparable 
to those of powerful ERI and ERII sources. There is no correlation 
between luminosity and spectral index. 

The redshift distribution for sources in the sample is compared 
to the distribution of the catalogue; the samples have similar dis¬ 
tributions, although the sources with z < 0.2 which are predicted 
to be present by the simulation are missing from the IOC sample. 
These low-redshift sources in the simulated sample are starform¬ 
ing sources, so the fact that they are missing from the sample is 
consistent with the finding that there are essentially no starforming 
sources in the IOC sample. 

The proportion of radio-loud sources in the IOC 15.7-GHz se¬ 
lected sample (> 94 percent) is significantly higher than the propor¬ 
tion in the ATESP 5-GHz selected sample (~ 60 percent), which has 
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a comparable flux density range. High frequency surveys are there¬ 
fore a very effective method of selecting sub-mjy radio-loud AGN, 
as the steep-spectrum starforming galaxies found in samples selec¬ 
ted at lower frequencies are not present. The radio galaxies in this 
sample will be discussed further in a following paper, where they 
will be split into high-excitation and low-excitation radio galaxies 
and the properties of these two classes of sources are compared. 
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APPENDIX A: MULTI-WAVELENGTH DATA 

Table A1 shows the multi-wavelength data and redshift estimates 
available for the 96 IOC sources studied in this paper. Table A2 
lists the origins of the spectroscopic redshifts used. 
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Table Al. Multi-wavelength properties of the 96 IOC sources studied in this paper. Optical and infrared magnitudes or flux densities are listed, along with redshift values. Redshifts from the Le Phare photometric 
redshift fitting described in Section 4 are given, as well as values from the Rowan-Robinson et al. (2013) and Fotopoulou et al. (2012) catalogues. The ‘best’ redshift value, used in the analysis, is also listed. All 
magnitudes are AB magnitudes. 
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5 = compact, match; 6 = compact, no match; 7 = IOC position only. 

b) Final redshift value, description is given in Section 4.3. 

c) Origin of final redshift value. 1 = spectroscopic (see Table A2) 2 = RR13, 3 = F12, 4 = Le Phare. 
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Table Al - continued 
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Table A2. References for the spectroscopic redshifts used in this work. 


AMI ID 

Spec, z 

Reference 

IOC J104320+585621 

0.35 

NED" 

IOC J104428+591540 

0.36 

NED 

IOC J104528+591328 

2.31 

NED 

IOC J104539+585730 

0.39 

NED 

IOC J104630+582748 

0.12 

CfA HectoSpec Spitzer Follow-Up by Huang, Rigopoulou et al. (in prep) 

IOC J104633+585816 

0.85 

NED 

IOC J104719+5821I4 

1.22 

NED 

IOC J104918+582801 

2.30 

CfA HectoSpec Spitzer Follow-Up by Huang, Rigopoulou et al. (in prep) 

IOC J104939+583530 

0.97 

SDSS-DR12, Alam et al. (2015) 

IOC J104943+571739 

0.59 

CfA HectoSpec Spitzer Follow-Up by Huang, Rigopoulou et al. (in prep) 

IOC J104954+570456 

0.53 

CfA HectoSpec Spitzer Follow-Up by Huang, Rigopoulou et al. (in prep) 

IOC J105028+574522 

0.07 

NED 

IOC J105039+572339 

1.44 

NED 

IOC J105039+585118 

0.37 

SDSS-DR12'', Alam et al. (2015) 

IOC J105104+5754I5 

1.67 

SDSS-DR12, Alam et al. (2015) 

IOC J105128+570901 

0.54 

SDSS-DR12, Alam et al. (2015) 

IOC J105132+571114 

0.32 

NED 

IOC J105148+573245 

0.99 

NED 

IOC J105206+574111 

0.46 

NED 

IOC J105225+573323 

0.61 

SDSS-DR12, Alam et al. (2015) 

IOC J105425+573700 

0.32 

ITP 2010 HerMES Follow-Up by Perez-Fournon, Page et al. (in prep.) 

IOC J105527+571607 

0.49 

SDSS-DR12, Alam et al. (2015) 

IOC J105550+570407 

0.49 

SDSS-DR12, Alam et al. (2015) 

IOC J105653+580342 

0.60 

SDSS-DR12, Alam et al. (2015) 


Notes: 

a) NASA/IPAC Extragalactic Database (extracted June 2013). 

b) Sloan Digital Sky Survey Data Release 12. 
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